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; ABSTRACT 

\^ , We present high-precision photometry of three transits of the extrasolar planetary system 

■ WASP-2, obtained by defocussing the telescope, and achieving point-to-point scatters of be- 

\ tween 0.42 and 0.73 mmag. These data are modelled using the JKTEBOP code, and taking 

into account the light from the recently-discovered faint star close to the system. The physical 
properties of the WASP-2 system are derived using tabulated predictions from five different 
sets of stellar evolutionary models, allowing both statistical and systematic errorbars to be 
. specified. We find the mass and radius of the planet to be Afb = 0.847 ± 0.038 ± 0.024 Mjup 

^ ; and i?b = 1.044 ± 0.029 ± 0.015 Rjup. It has a low equilibrium temperature of 1280 ± 21 K, 

in agreement with a recent finding that it does not have an atmospheric temperature inversion. 
The first of our transit datasets has a scatter of only 0.42 mmag with respect to the best-fitting 
light curve model, which to our knowledge is a record for ground-based observations of a 
transiting extrasolar planet. 

Key words: stars: planetary systems — stars: individual: WASP-2 — stars: binaries: eclipsing 



1 INTRODUCTION 

Aside from our own Solar system, currently the only way to mea- 
sure the mass and radius of a planet is to study one which tran- 

* Based on data collected by MiNDSTEp with the Danish 1 .54 m telescope 

at the ESO La Silla Observatory 

t Royal Society University Research Fellow 



sits its parent star. Almost all known transiting extrasolar planets 
(TEPs) are gas giants with radii similar to Jupiter's but who or- 
bit much closer to their parent stars. The masses and radii, and 
thus surface gravities and mean densities, of these objects can be 
determined to high precision, yielding in turn constraints on their 
chem i cal composition an d internal structure ( Bodenheim er et al] 
I2OO3I : iBaraffe et al.ll2008l : iFortnev et alj|2007l boosi) . The flux ra- 
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Table 1. Log of the observations presented in this work. N^^^ is the number of observations and 'Moon ilium.' is the fractional illumination of the Moon at 
the midpoint of the transit. 



Transit 


Date 


Start time 


End time 




Exposure 


Filter 


Airmass 


Moon 


Aperture 


Scatter 






(UT) 


(UT) 




time (s) 






ilium. 


sizes (px) 


(mmag) 


1 


2009 06 03 


05:44 


10:35 


123 


90-100 


Rc 


1.66 -> 1.23 ~s> 1.48 


0.816 


26, 33, 100 


0.422 


2 


2009 08 1 1 


02:27 


05:06 


97 


50.0 


Rc 


1.34 1.23 1.30 


0.761 


23, 40, 70 




3 


2009 08 24 


01:08 


05:07 


83 


120.0 


Rc 


1.42 1.23 1.45 


0.179 


30, 35, 85 


0.473 


4 


2009 11 11 


17:00 


20:11 


63 


130.0 


Rc 


1.28^2.17 


0.273 


20, 30, 45 


0.730 



tio of a planet to its parent star can also be obtained through mea- 
surements of the depth of its occultation (secondary eclipse) at a 
range of wavelengths, delivering constraints on the st ructure of ex- 
tended and highly irra diated atmospheres (e.g. iKnutson et alj2008l : 
lAnderson et"al]|2010l) . 

To measure the physical properties of a TEP we need a good 
light curve covering one or more transits, radial velocity (or astro- 
metric) observations of its parent star at a range of orbital phases, 
and one additional constraint usually supplied by forcing the prop- 
erties of the parent star to match predictions from theoretical stellar 
models. Detailed error budgets of the analysis process ( Southworth 
I2OO9I) show that the dominant contributor to the uncertainties is 
almost always the light curve. We have therefore embarked on a 
project to obtain high-quality light curves of transiting systems and 
thus measure their physical properties as accurately and homoge- 
neously as possible. 

Our approach to obtaining transit light curves is to heavily de- 
focus a telescope, spreading the light from each star over thousands 
of pixels on the CCD, and take images with comparatively long ex- 
posure times of up to 120 s. Because the light from each star is 
distributed over thousands of pixels, it is possible to achieve both a 
very low Poisson noise and minimal flat-fielding errors. Using long 
exposure times mean we spend less time overall reading out the 
CCD, so are able to detect more light from the star and decrease 
the noise due to scintillation. Defocussing the telescope means our 
observations are subject to a higher background light level, but this 
is not important for the bright stars around which most known TEPs 
have been found. In this work we present defocussed observations 
of the planetary system WASP-2, obtained using two telescopes and 
achieving scatters ranging from 0.42 to 0.73 mmag per datapoint. 

The transiting nature of WASP-2 was discovered by 
Collier Cam eron et al. ( 2007b) (see also ICollier Cameron et al.l 



2007af) in photome try from the SuperWASP-North telescope 
jPoUacco et ai]|2006l) . The planet WASP-2 b is slightly larger and 
less massive than Jupiter, whilst the star WASP-2 A is rather 
smaller and cooler than the Sun. Follow-up pho tometry of one tran- 
sit was presen ted and analysed by Cha rbonneau et all j2007l). an d 
reanalysed by [Torres et all ( |2008() and ISouthworthI bOOSi, |2009|) . 



Another transit has also been observe d with the W i lliam Her- 
schel Telescope jHrudkova l2009h . iTriaud et all ( 1201 Oh has 
presented a study of WASP-2 which incorporates spectroscopic 
observations of the Rossiter-McLaughlin effect, spectroscopic esti- 
mates of the effective temperature and metal ab undances of the host 
star, a nd a re-analysis of th e light curves from Charbonneau et al] 
( l2007t) and lHrudkovaetal.1 11)09). Whe atlev et alj ( I2OI0I) has re- 
cently presented infrared Spitzer observations of the secondary 
eclipse of the system. 

The study of WASP-2 presents an additional complication 
in the recently-disco vered presence o f a faint companion star to 
the WASP-2 system. iDaemgen et alj ( l2009i) used the high-speed 



"lucky-imaging" camera AstraLux to obtain images of the field of 
WASP-2 in the SDSS i and z passbands and with a high spatial 
resolution. These authors discovered a faint point source within 
0.76" of WASP-2, and measured magnitude differences of Ai — 
4.095 ± 0.025 and Az = 3.626 ± 0.016 mag with respect to our 
target star. These correspond to a light ratio of 0.0354 ± 0.0005 
and 0.0230 ± 0.0005 in the two filters, respectively. Because this 
faint star was within the point spread function (PSF) of all previous 
photometric observations, studies which relied on these data will 
be systematically incorrect by a small amount. In Section [T2l we 
outline our method for accounting for the light from this faint star. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed three transits of WASP-2 in 2009 June and August, 
using the 1 .54 m Danish Telescop43 at ESO La Silla, Chile. We also 
observed one transit on the night of 2009 November 1 1 using the 
1.52 m G. D. Cassini Telescopqjat Loiano Observatory, Italy. The 
second of our transit observations was strongly affected by cloud, 
so we do not include these data in our analysis. The fourth transit 
(Cassini) was also taken in non-photometric conditions, causing the 
observations to be noisier than intended. An observing log is given 
in Tabled 

Observations taken with the Danish telescope used the 
DFOSC focal-reducing imager, which has a plate scale of 
0.39" px^^. The CCD was not binned, but was windowed down 
by varying amounts to decrease the readout time. All observations 
were done through a Cousins R filter, and the amount of defo- 
cussing was adjusted until the peak counts per pixel from WASP- 
2 were roughly 25 000 above the sky background, resulting in 
doughnut- shaped PSFs with diameters ranging from 30 to 45 pix- 
els. The pointing of the telescope was maintained using autogu- 
iding, and we did not change the amount of defocussing during 
individual observing sequences. The observing sequence from the 
Cassini telescope used the BFOSC focal-reducing imager (plate 
scale 0.58" px~^) and the same approach as for the Danish tele- 
scope. 

Several images were taken with the telescopes properly fo- 
cussed, to verify that there were no faint stars within the defocussed 
PSF of WASP-2. The closest detectable star is 4.7 mag fainter and 
at a distance of 18", and its PSF did not overlap with our target star. 
We conclude that no stars interfere with the PSF of WASP-2 aside 
from the one at 0.76" discovered by Daemgen et al. (2009). 

Data reduction was performed as in previous papers 



1 For information on the 1.54 m Danish Telescope and DFOSC see: 



Ihttp : / /www . e so. org/ sci/ facilities /las ilia/ teles cope s/dlp5 / 

^ Information on the 1.52m Cassini Telescope and BFOSC can be found 
at,http : / / www .bo.astro.it/loiano/ 
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Figure 1. Light curves of WASP-2 from the three nights of observations 
included in our analysis. The en'or bars have been scaled to give x'^,^ = 1.0 
for each night, and some are smaller than the symbol sizes. 



Table 2. Exceipts of the Hght curve of WASP-2. The full dataset will be 
made available as a supplement to this paper and at the CDS. 



Telescope 


HJD 


Diff. mag. 


Uncertainty 


Danish 


2454985.742452 


0.0009568 


0.0004018 


Danish 


2454985.944273 


-0.0003249 


0.0004759 


Danish 


2455054.609274 


0.0025309 


0.0005716 


Danish 


2455054.718127 


0.0246546 


0.0016991 


Danish 


2455067.553072 


-0.0011116 


0.0003926 


Danish 


2455067.719132 


-0.0001831 


0.0003860 


Cassini 


2455147.207773 


0.0012495 


0.0006328 


Cassini 


2455147.341286 


-0.0009452 


0.0005788 



dSouthworth et al.ll2009allbll3) . using a custom pipeline written in 
the IDl ^pro^ranm iing language and accessing the DAOPHOT pack- 
age (Stetson 1987) to perform aperture photometry with the APER 
□ routine. The apertures were placed by eye. For transits 1 and 3 



their positions were held fixed throughout the observing sequence. 
For transits 2 and 4 the pointing was monitored by cross-correlating 
each image with a reference image, and the apertures shifted to fol- 
low the stars. We tried a wide range of aperture sizes and retained 
those which gave photometry with the lowest scatter; in each case 
we find that the shape of the light curve is very insensitive to the 
aperture sizes. 

Up to 11 comparison stars were measured on each image, and 
checked for short-period variability. We then selected the five best 
and combined them into a weighted ensemble which minimised the 
scatter in observations before and after transit. Simultaneously with 
this procedure we fitted a straight line to the data outside transit to 
normalise the light curve and remove slow magnitude drifts due to 
atmospheric effects. 

We also tried flat-fielding and debiassing the observations be- 
fore measuring aperture photometry. We found that our treatment 
of the CCD bias level had a negligible effect, and that flat-fielding 
the data actually made things worse. We therefore debiassed the 
data but did not flat-field them. The final light curves are shown in 
Fig.[T]and tabulated in Table[2] The scatter in the final light curves 
varies from 0.42 to 0.73 mmag. To our knowledge the first tran- 
sit dataset represents the most precise light curve of a point source 
ever obtained from a ground-based telescope. 



3 LIGHT CURVE ANALYSIS 

The analysis of our l ight curves was performed usin g the approach 
discussed in detail bv lSouthworti] i2008ll2009l.l2010h . so our results 
are consistent with the homogeneous analysis of TEPs presented in 
those works. I n short, the li ght curves were modelled using the JK- 
TEB0i|3 code dSouthworth et al. 2004a b), which is based on the 
EBOP program jPopper & Etzell98 1. ; ,Etzellll98ll : iNelson & David 
I1972I) where the two components of a binary system are simulated 
using biaxial spheroids. The asphericity of the components is gov- 
erned by the mass ratio, for which we used the value 0.0009. Large 
changes in this value have a negligible effect on our results, partly 
because in a system with a circular orbit the times of eclipse are 
when the sky projections of the components suffer the smallest dis- 
tortion. 



3.1 Period determination 

As a first step in the analysis of WASP-2 we have performed a 
fit to each of our three light curves separately using JKTEBOP. The 
measurement errors of our photometry (which come from the APER 
routine) were found to be underestimated so were rescaled to give a 
reduced y'^ of x'j = 1.0 with respect to the best-fitting model (see 
iBruntt et al.ll2006h . 

These individual fits also yielded one time of minimum 
light each. The uncertainties in these values were estimated us- 
ing Monte Carlo simulations. We have augmented them with thre e 
times of minimum light ta ken from Collier Cameron et al. Il l2007bl) . 
ICharbonneau et alj [2007h and lHrudkova et al.l ilOO^ . and with an 
additional 22 minimum timings obtained by amateur astronomers 
and made available on the AXJ^ and TRESC^] websites. A small 



The acronym IDL stands for Interactive Data Language and is a 
trademark of ITT Visual Information Solutions. For further details see 
[http : / /www . ittvis . com/Product Services /IDL . aspx 
* APER is part of the ASTROLIB subroutine library distributed by NASA. 
For further details see http : / / idlastro . gsf c . nasa . gov/. 



^ JKTEBOP is written in FORTRAN77 and the source code is available at 
http : / / www . astro. keele. a c ■ uk/~ jkt/| 

^ Amateur Exoplanet Archive, http ; / /brucegary . net /AXA/ x . html 
^ The TRansiting ExoplanetS and CAndidates (TRESCA) website can be 
found at, http://var2. astro. cz/EN/tresca/ index . php 
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Figure 2. Plot of the residuals of the timings of mid-transit of WASP-2 versus a linear ephemeris. Some error bars are smaller than the symbol sizes. The 
timings obtained by amateur astronomers are plotted using open circles. 



Table 3. Times of minimum light of WASP-2 and their residuals versus the 

ephemeris derive d in this vyork. 

Refere nce s: (1) ICoUier Cameron et al J j2007bl) ; (2) ICharbonneau et alj 
)2007l) : (3) lHrudkova et al.l )2009h ^ 



Time of minimum 


Cycle 


Residual 


Reference 


(HJD - 2400000) 


no. 


(HJD) 




53991.5146 ±0.0044 


0.0 


0.00005 


1 


54008.73205 ± 0.00028 


8.0 


-0.00027 


2 


54357.39254 ±0.00016 


170.0 


0.00034 


3 


54985.84071 ±0.00009 


462.0 


-0.00015 


This work 


55067.62554 ±0.00011 


500.0 


0.00027 


This work 


55147.25721 ±0.00017 


537.0 


-0.00025 


This work 


54260.5420 ±0.0010 


125.0 


-0.00023 


Gregorio (AXA) 


54277.7580 ±0.0020 


133.0 


-0.00200 


Gary (AXA) 


54279.9110 ±0.0010 


134.0 


-0.00122 


Foote (AXA) 


54316.4991 ±0.0015 


151.0 


-0.00089 


Vanmunster (AXA) 


54320.8053 ±0.0010 


153.0 


0.00087 


Foote (AXA) 


54333.7186 ± 0.0020 


159.0 


0.00084 


Sheridan (AXA) 


54348.7760 ±0.0020 


166.0 


-0.00731 


Sheridan (AXA) 


54357.3906 ±0.0010 


170.0 


-0.00160 


Poddany (AXA) 


54357.3914 ±0.0010 


170.0 


-0.00080 


Ohlert (AXA) 


54710.35882 ±0.00100 


334.0 


0.00231 


Kocian (TRESCA) 


54757.7052 ±0.0006 


356.0 


-0.00018 


Gary (AXA) 


54953.55627 ± 0.00083 


447.0 


-0.00126 


Brat (TRESCA) 


54981.53709 ±0.00044 


460.0 


0.00068 


Sauer (TRESCA) 


55009.5168 ±0.0018 


473.0 


0.00151 


Gonzalez (AXA) 


55024.5839 ±0.0020 


480.0 


0.00306 


Lobao (AXA) 


55026.7327 ±0.0006 


481.0 


-0.00036 


Norby (AXA) 


55037.4900 ±0.0014 


486.0 


-0.00417 


Salom (AXA) 


55037.4945 ±0.0010 


486.0 


0.00033 


Naves (AXA) 


55065.4735 ±0.0008 


499.0 


0.00045 


Gregorio (AXA) 


55069.7799 ±0.0015 


501.0 


0.00241 


Garlitz (AXA) 


55074.08138 ±0.00050 


503.0 


-0.00055 


Cao (TRESCA) 


55095.6050 ±0.0009 


513.0 


0.00085 


Cordiale (AXA) 



number of AXA and TRESCA timings were rejected due to noisy 
data or incomplete coverage of the transit. 

Armed with these times of minimum light we have redeter- 
mined the orbital ephemeris of the WASP-2 system to be 

To = HJD 2 453 991.51455(17) + 2.15222144(39) x E 

where E is the number of orbital cycles after the reference epoch 



and quantities in parentheses denote the uncertainty in the final 
digit of the preceding number. The of the straight line fitted 
to the timings is equal to 1.53, so the ephemeris uncertainties were 
multiplied by the square-root of this number to yield the values 
quoted above. A plot of the residuals around the fit is shown in 
Fig.|2] Including the 22 amateur epochs roughly doubles the preci- 
sion of the ephemeris. 

There were several cases during the 2009 MiNDSTEp observ- 
ing season, which is primarily oriented to studying microlensing 
events towards the Galactic bulge, where photometric observations 
appeared to occur slightly earlier than expected. The timestamps 
iinprinted in the DFOSC fits file headers are currently under inves- 
tigation. The Julian Dates for the DFOSC observations presented 
in the current work should therefore be treated with caution, as 
should the observations of WASP- 18 studied bv lSouthworth et al] 
l l2009ch . One of the reasons for obtaining a further transit obser- 
vation of WASP-2 with the Cassini telescope was to improve its 
orbital ephemeris and allow an external check on the DFOSC tim- 
ings. We find that the two DFOSC transits occur 12 ± 8 s later than 
would be expected, implying that there is no problem with their 
timestamps. We have therefore included these timings in our pe- 
riod study. 



3.2 Accounting for the light from the faint star 

Light from the faint star, discovered by iDaemgen et al.l ( l2009h 
0.76" away from WASP-2, must be accounted for to avoid sys- 
tematically underestimating the radius of the planet. The PSF of 
WASP-2 in our observations includes essentially all the light from 
the faint star, and even our focussed observations were unable to 
pick this star up due t o its fa intness and closeness to our target star. 

i Daemgen etal] ( l2009h obtained light ratios between the 
WASP-2 system and the faint star in the SDSS i and z bands. 
We have extrapolated these light ratios to the Rc band using 
synthetic spec t ra cal culated with the ATLAS9 model atmospheres 
l lKurucz|[l979ll 19931) and filter response functions taken from the 
Isaac Newton Group websit^fl We find an effective temperature 
for the faint star of T^s = 3780 ± 46K and a light ratio of 



Isaac Newton Group: www . ing .iac.es 
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Figure 3. Phased light curves of WASP-2 from the Da nish telescope (top), 
the Ca ssini telescope (middle grey points) and from ICharbonneau et al] 
(lower), compared to the best fits found using JKTEBOP and the 
quadratic LD law. The residuals of the fits are plotted at the bottom of the 
figure, offset from zero. 

l-R. = 0.0176 ± 0.0007, where the errorbars include contribu- 
tions from uncertainties in the i and z light ratios and in the Tcff 
of WASP-2 A. There will be a systematic error in from the use 
of theoretical model atmosphere calculations, but this will be very 
small as the effective wavelength of Rc is close to that of i. The 
systematic error in the T^s of the faint star will be substantially 
larger, but can be ignored because we do not use this Tcff measure- 
ment in our analysis. 

The light from the faint star can be straightforwardly included 
in our photometric analysis, as JKTEBOP (and other light curve 
codes) allow for thi s possibility under the moniker 'third light' 
(symbolised by L^). ISouthworthI ( I2OI0I) modified JKTEBOP to al- 
low L3 to be included as a fitted parameter constrained by a mea- 
sured value. For our defocussed observations, which were all taken 
through t he Rc filter, we adopt L3 = 0.0176 ±0.0007. The obser- 
vations of ICharbonneau et alj ( l2007h were obtained using an SDSS 
z filter, fo r which a magn i tude d ifference (light ratio) was given 
directly bv lOaemgen et al.l ( |2009|) . 



3.3 Light curve modelling 

Our light curves have been modelled using the JKTEBOP code. The 
main parameters of the model are the fractional radii of the star and 
planet, rx and rb (the radii of the components divided by the semi- 
major axis), and the orbital inclination, i. va and rb are actually 



incorporated as their sum, rx+rh, and ratio, k, as these parameters 
are more weakly correlated. 

The two Danish telescope datasets were solved simultane- 
ously and the Cassini telescope data were studied separately. 
We have also rev i sited the z-band light curve presented by 
Charbonne au et al] ( l2007h . and this solution supersedes that of 
Southworth ( 2008) which was presented before the faint star nearby 
WASP-2 was discovered. 

Limb darkening (L D) was included using any of five LD laws 
(see ISouthworthI I2OO8I) and with theoretical coeffici ents (ma and 
va) ob t ained by bilinear inte rpol ation in the tables of IVan Hammd 
( ll993h . lClaretl j2000ll2004al) and lClaret & Hauschildtl ilOOA So- 
lutions were obtained with ua and va both fixed, both included as 
fitted parameters, and with it a fitted and va perturbed by ±0.10 
on a flat distribution. For our final solutions we adopt the last of 
these three possibilities, as the data cannot support the extraction 
of both LD coefficients. This does not cause a significant depen- 
dence on stellar th eory because the two LD coefficients are very 
strongly correlated jSouthworth et alj2007ah . 

We assessed the uncertainties of the fitted parameter s using 
1000 Monte C arlo simulations for each solution ( Southw orth et al. 
2004cl l2005h . The residual-permutation method jjenkins et al.l 
2002h was also used to account for correlated noise, which was 
found to be unimportant for all sets of observations. The results for 
each light curve are given in Tables|4][5]and[6] 

The final value quoted for each light curve parameter is the 
unweighted mean of the values from the four nonlinear LD law so- 
lutions. Its uncertainty is the larger of two possibilities: the largest 
of the individual Monte Carlo uncertainties with the largest devi- 
ation from the mean value added in quadrature; or the uncertainty 
found by the residual-permutation method. The best fits are shown 
in Fig.[3]and their parameters are tabulated in Table|7] The agree- 
ment is unusually good, with xt values of between 0. 1 and 0.3 in 
the final weighted-mean values. Table[8] contains literature values 
of the light curve parameters; a direct comparison with our results 
in Table|7] shows not only a good agreement but also a substantial 
decrease in the errorbars. A quick check of the Monte Carlo sim- 
ulation output shows that L3 is almost uncorrelated with the other 
fitting parameters. The inclusion of third light in the light curve so- 
lutions does not cause any difficulties in this case because its value 
is small and well-determined. 



4 THE PHYSICAL PROPERTIES OF WASP-2 

In order to determine the physical properties we need the light 
curve parameters determined above, the velocity amplit ude of the 
star (i^A = 153.6 ± 3.0 km s~^ ; iTriaud et al] |2010|) . and one 
additional constraint. We supply this extra constraint by forcing 
the physical properties of the star to match predictions from the- 
oretical stellar models, guided by the observed effective temper- 
ature (Tcff = 5150 ± 80 K; Jriaud et al. 2010) and metallicity 
([§] = -0.08 ± 0.08; ITriaud etalj|201 0). For the parameter 
which governs the solution process we use the velocity amplitude 
of the planet, Kh- 

For the theoretical models we use five differe nt sets of tabula- 
tions: Claret (jcfaret 2004b, 2005, 2006, 2007J, Y '^ llDemaraue et al. 
200?), Teramo (Pietrinferni et a l.ll2004l) . VRSS jVandenBerg et al] 
I200Q) and DSEP tDotter et al.ii2008l). A detailed des cription of this 
solution process can be found in ]Southworthl ( |2009|) and a discus- 
sion of the five stellar mo dels and the altern ative empirical mass- 
radius relation is given by ISouthworthI ( [ioToi) . The use of five dif- 
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Table 4. Parameters of the JKTEBOP best fits of the Danish telescope light curve, using different approaches to LD. For each part of the table the upper 
quantities are fitted parameters and the lower quantities are derived parameters. Tq is given as HJD — 2454000.0. The light curve contains 211 datapoints. 



Linear LD law Quadratic LD law Square-root LD law Logarithmic LD law Cubic LD law 



All LD coefficients fixed: 



rA + rb 0.1391±0.0016 0.1387 ± 0.0016 0.1392±0.0016 0.1394 ± 0.0016 0.1424±0.0014 

k 0.13342 ±0.00054 0.13217 ± 0.00058 0.13269 ±0.00056 0.13249 ±0.00056 0.1 3338 ±0.00041 

j(deg.) 84.880 ±0.1 12 84.960 ±0.115 84.902 ±0.107 84.886 ±0.106 84.568 ±0.096 

«A 0.60 fixed 0.50 fixed 0.35 fixed 0.72 fixed 0.40 fixed 

UA 0.20 fixed 0.42 fixed 0.21 fixed 0.10 fixed 
To 985.840869 ±0.000071 985.840868 ± 0.000072 985.840868 ± 0.000073 985.840868 ± 0.000071 985.840869 ±0.000067 

rA 0.1227±0.0014 0.1225 ± 0.0014 0.1229 ± 0.0013 0.1231 ± 0.0014 0.1256 ±0.0012 

rb 0.01637 ±0.00024 0.01619 ± 0.00025 0.01631 ±0.00023 0.01631 ± 0.00024 0.01676 ±0.00020 

(T(mmag) 0.4833 0.4857 0.4843 0.4843 0.4919 

xl 1.1452 1.1580 1.1505 1.1505 1.1874 



Fitting for the linear LD coefficient and perturbing the nonlinear LD coefficient: 

rA + rb 0.1393±0.0019 0.1397 ± 0.0020 0.1398 ± 0.0020 0.1398 ±0.0020 0.1398±0.0020 

k 0.1 3350 ±0.00066 0.13259 ± 0.00068 0.13293 ±0.00064 0.13262 ± 0.00071 0.1 3299 ±0.00071 

j(deg.) 84.85 ±0.15 84.84 ±0.16 84.83 ±0.16 84.84 ±0.16 84.83 ±0.16 

itA 0.585 ±0.043 0.437 ±0.070 0.312 ±0.061 0.692 ±0.062 0.534 ±0.055 

v\ 0.20 perturbed 0.42 perturbed 0.21 perturbed 0.10 perturbed 

To 985.840869 ±0.000067 985.840868 ± 0.000071 985.840869 ± 0.000070 985.840869 ± 0.000073 985.840869 ± 0.000071 

rA 0.1229±0.0017 0.1234 ± 0.0017 0.1234± 0.0017 0.1234 ± 0.0017 0.1234 ± 0.0018 

rb 0.01641 ±0.00029 0.01636 ± 0.00029 0.01641 ±0.00030 0.01637 ± 0.00030 0.01640 ±0.00030 

cr(mmag) 0.4831 0.4839 0.4837 0.4840 0.4837 

xl 1-1499 1.1539 1.1529 1.1541 1.1528 



Table 5. As Table|4]but for the Cassini data. To is given as HJD — 2455000.0 and the light curve contains 63 datapoints. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


All LD coefficients fixed: 










rA + rb 


0.1382± 0.0041 


0.1371 ±0.0046 


0.1379 ±0.0042 


0.1380 ±0.0044 


0.1415 ±0.0035 


k 


0.1320 ±0.0015 


0.1305 ±0.0017 


0.1311 ±0.0016 


0.1309 ±0.0016 


0.1324±0.0011 


i (deg.) 


85.01 ±0.29 


85.15 ±0.32 


85.07 ±0.30 


85.05 ±0.31 


84.69 ±0.23 


MA 


0.60 fixed 


0.50 fixed 


0.35 fixed 


0.72 fixed 


0.40 fixed 






0.20 fixed 


0.42 fixed 


0.21 fixed 


0.10 fixed 


To 


147.25724 ±0.00017 


147.25725 ±0.00018 


147.25724 ±0.00017 


147.25724 ±0.00018 


147.25721 ±0.00017 


rA 


0.1221 ±0.0035 


0.1212 ±0.0038 


0.1219 ±0.0036 


0.1221 ±0.0038 


0.1249 ±0.0030 


rb 


0.01611 ±0.00062 


0.01582 ± 0.00070 


0.01598 ±0.00063 


0.01598 ±0.00066 


0.01654 ±0.00050 


a (mmag) 


0.7385 


0.7417 


0.7387 


0.7371 


0.7260 




1.0248 


1.0344 


1.0255 


1.0206 


0.9875 


Fitting for the linear LD coefficient and perturbing the nonlinear LD coefficient: 








0.1426 ±0.0049 


0.1427 ±0.0052 


0.1427 ±0.0051 


0.1427 ±0.0049 


0.1427 ±0.0050 


k 


0.1330±0.0014 


0.1321 ±0.0014 


0.1324 ±0.0014 


0.1322 ±0.0013 


0.1325 ±0.0014 


i (deg.) 


84.58 ±0.41 


84.58 ± 0.42 


84.58 ±0.41 


84.58 ±0.38 


84.58 ±0.41 


UA 


0.39±0.18 


0.23 ±0.21 


0.12±0.19 


0.49 ±0.19 


0.34 ±0.20 


VA 




0.20 perturbed 


0.42 perturbed 


0.21 perturbed 


0.10 perturbed 


To 


147.25721 ±0.00018 


147.25721 ±0.00017 


147.25721 ±0.00017 


147.25721 ±0.00017 


147.25721 ±0.00017 


rA 


0.1258 ±0.0042 


0.1261 ±0.0045 


0.1260 ±0.0044 


0.1260 ±0.0042 


0.1260 ±0.0043 


rb 


0.01673 ±0.00063 


0.01666 ±0.00068 


0.01668 ±0.00065 


0.01666 ±0.00065 


0.01669 ±0.00064 


a (mmag) 


0.7251 


0.7245 


0.7246 


0.7245 


0.7246 




1.0021 


1.0003 


1.0007 


1.0002 


1.0007 



ferent model sets allows the assessment of their interagreement 
and thus the systematic error engendered by relying on theoreti- 
cal predictions. In reality this will yield only a lower limit on the 
systematic errors, as the different models have a lot of common- 
ality in their input physics. There are two parameters which do 
not depend on the stellar model tabulations - the surface grav- 
ity and equilibrium temperature of the planet jSouthworth et al.l 



l2007bl : ISouthworthll201(]|) - a nd one for which the dependenc e is 
negligible - the stellar density dSeager & Mallen-Omela3l2003h . 

The results of the above analysis are estimates of the physical 
properties of the WASP-2 system, given in Table|9]for each stellar 
mo del set and for the e mpirical stellar mass-radius relation defined 
bv lSouthworthl l2008). The mass, radius, surface gravity and den- 
sity of the star are denoted by A/a, Ra, loggA and pA, and of 
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Linear LD law Quadratic LD law Square-root LD law Logarithmic LD law Cubic LD law 



All LD coefficients fixed: 





0.1 390 ± 0.0049 


0.1397 ±0.0045 


0.1 399 ±0.0044 


0.1399 ±0.0047 


0.1416 ±0.0041 


k 


0.1343 ±0.0016 


0.1331 ±0.0014 


0.1337±0.0013 


0.1334 ±0.0014 


0.1340 ±0.0011 


i (deg.) 


85.01 ±0.35 


84.97 ±0.31 


84.92 ±0.31 


84.93 ±0.33 


84.73 ±0.28 


UA 


0.54 fixed 


0.29 fixed 


0.13 fixed 


0.61 fixed 


0.30 fixed 


VA 




0.30 fixed 


0.53 fixed 


0.26 fixed 


0.10 fixed 


To 


8.73204 ±0.00021 


8.73204 ± 0.00022 


8.73205 ± 0.00020 


8.73205 ± 0.00022 


8.73205 ± 0.00020 


rA 


0.1225 ±0.0042 


0.1233 ±0.0039 


0.1234 ±0.0038 






rh 


0.01646 ± 0.00072 


0.01641 ±0.00064 


0.01649 ±0.00064 


0.01646 ±0.00067 


0.01674 ±0.00057 


u (mmag) 


1.8878 


1.8873 


1.8869 


1.8869 


1.8891 




0.9975 


0.9970 


0.9965 


0.9966 


0.9989 


Fitting for the linear LD coefficient and perturbing the nonlinear LD coefficient: 






rA + r-b 


0.1395 ±0.0057 


0.1400 ±0.0060 


0.1399 ±0.0055 


0.1400 ±0.0060 


0.1 397 ±0.0058 


k 


0.1344 ±0.0016 


0.1332 ±0.0014 


0.1337 ±0.0015 


0.1335 ±0.0015 


0.1339 ±0.0016 


i (deg.) 


84.93 ± 0.46 


84.91 ±0.46 


84.92 ± 0.45 


84.91 ±0.47 


84.93 ±0.47 


UA 


0.47 ±0.15 


0.24 ±0.17 


0.13±0.16 


0.60±0.16 


0.42 ±0.16 


VA 




0.30 perturbed 


0.53 perturbed 


0.26 perturbed 


0.10 perturbed 


To 


8.73204 ± 0.00022 


8.73205 ± 0.00022 


8.73205 ±0.00021 


8.73205 ±0.00022 


8.73205 ±0.00021 


TA 


0.1229 ± 0.0049 


0.1236±0.0051 


0.1234 ±0.0048 


0.1236 ±0.0052 


0.1232 ±0.0050 


rh 


0.01652 ±0.00081 


0.01646 ±0.00078 


0.01650 ±0.00077 


0.01649 ±0.00083 


0.01650 ±0.00081 


a (mmag) 


1.8869 


1.8869 


1.8869 


1.8869 


1.8869 




0.9988 


0.9989 


0.9989 


0.9989 


0.9988 



Table 7. Parameters of the fits to the three light curves of WASP-2. The final parameters are the weighted mean of the result for the three light curves. 





This work (Danish) 


This work (Cassini) 


This work (Charbonneau) 


This work (final) 


''A + fb 


0.1398 ±0.0024 


0.1427 ±0.0057 


0.1399 ±0.0067 


0.1402 ±0.0021 


k 


0.1328 ±0.0008 


0.1323 ±0.0015 


0.1336 ±0.0017 


0.1328 ±0.0007 


in 


84.84 ±0.19 


84.58 ±0.47 


84.92 ± 0.52 


84.81 ±0.17 


rA 


0.1234±0.0021 


0.1260 ±0.0050 


0.1234 ±0.0058 


0.1237 ±0.0018 


rb 


0.01638 ±0.00036 


0.01667 ±0.00076 


0.01649 ±0.00092 


0.01644 ±0.00031 



These data were analysed with the JKTEBOP code, taking 
into account the contaminating light from a fainter companion 
star within 0.76" of the sky position of the WASP-2 system. 
With the addition of the observed effective temperature, metal- 
licity and velocity amplitude of the host star, plus several sets of 
theoretical stellar model predictions, we have obtained the physi- 
cal properties of the system. The planet WASP-2 b has a mass of 
0.847 ± 0.038 ± 0.024 Mjup (statistical and systematic errors) and 
a radius of 1.044 ± 0.029 ± 0.015 Rjup, which is typical for tran- 
siting gas giants. Its host star is a 0.8 M© and 0.8 R© Kl dwarf 
with a large and largely unconstrained age. The equilibrium tem- 
perature of the plane t, 1280 ± 21 K, puts it firmly into the pL class 
i lFortnev et al]|2008h . 

The high precision of our mass and radius measurements for 
WASP-2 b means this planet is now one of the best-understood 
TEPs. This is partly due to the small scatter in our light curves, but 
also to the physical configuration of the system. Transiting plan- 
ets with relatively low orbital inclinations suffer from less correla- 
tion between inclination and the stellar radius, resulting in better- 
defined light curve solutions. Low-inclination TEPs also spend a 
larger fraction of their time in the partial phases of their eclipses, 
making transit timings more accurate. We find uncertainties of as 
little as 8 s in the times of the three transits presented here: WASP-2 
would be an excellent candidate for a transit timing study, and such 



the planet by Afb, Rb, Qb and pb- The orbital semimajor axis (a), 
equilibrium temperature surrogate 




and lSafronovl (Il972l) number (0) are also determined. 

Tablel lOlcontains the final physical properties from this work, 
calculated as the unweighted mean of the individual results for the 
five stellar models. The statistical error for each parameter is the 
largest from the models, and the systematic error the standard devi- 
ation of the individual model results. TablellOlalso contains results 
from published studies of WASP-2, which are in good agreement 
with our results but are less precise. 



5 SUMMARY 

Obtaining a high-precision light curve is a crucial part of mea- 
suring the physical properties of transiting extrasolar planets. In 
this work we have presented photometric observations obtained 
through a heavily defocussed telescope and covering three transits 
of WASP-2. The first of these observing sequences has a scatter of 
only 0.42 mmag around the best-fitting model, which to our knowl- 
edge is a new record for a ground-based light curve of planetary 
transit. 
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Table 8. Light curve parameters of WASP-2 taken from the hterature, for co mpaiison to the results presented in Table|2] Quantities without uncertainties have 
been calculated from other parameters given in that study. Only the study bv lDaemgen et alj j2009l) accounts for the light from the faint companion star. 



Collier Cameron Charbonneau ToiTes et al. Southworth Daemgen+ et al. Triaud et al. Wheatley et al. 

et al. (2007b) et al. (2007) (2008) (2008) (2009) (2010) (2010) 

'■A+'-b 0.086 to 0.132 
k 0.119 to 0.140 

rA 0.123 
rb 0.0105 to 0.0202 



0.1309 ±0.0015 
84.74 ±0.39 



0.1423 

n 19^0+0.0031 

0.01648 



0.1409 ±0.0067 
0.1313 ±0.0017 

84.83 ±0.53 
0.1245 ±0.0058 
0.01635 ± 0.00093 



0.1420 
0.1332 ±0.0015 
84.80 ±0.39 
0.1250 
0.01701 



0.14155 

,+0.0010 
-0.0009 

84.73l°;;« 

i+0.0025 
-0.0024 
U.UlD/0_„ „„„4„ 



0.1342"* 



0.1248"* 



0.1334±0.0015 



Table 9. Estimates of the physical properties of WASP-2. In each case = 19.30 ± 0.82 m s-^, pA = 1-531 ± 0.067 pQ and T ' = 1280 ± 21 K. 



This work 
(Mass-radius) 



This work 
{Claret models) 



This work 
(y2 models) 



This work 
(Teramo models) 



This work 
(raS5 models) 



This work 
{DSEP models) 



Kb (kms-i) 
Ma(Mq) 

loggA (cgs) 
Mb ( Mjup) 

Rb ( Rjup) 
Pb (PJup) 

e 

a(AU) 
Age (Gyr) 



150.9 ±4.3 
0.777 ± 0.066 
0.798 ±0.028 
4.525 ±0.015 
0.828 ±0.050 
1.032 ±0.035 
0.754 ± 0.050 
0.0619 ±0.0024 
0.03000 ± 0.00085 



154.7 ±3.0 
0.838 ±0.048 
0.818 ±0.018 
4.536 ±0.016 
0.871 ±0.038 
1.058 ±0.029 
0.735 ± 0.046 
0.0604 ± 0.0020 
0.03076 ± 0.00060 



153.2 ±2.6 
0.813 ±0.041 
0.810±0.017 
4.531 ±0.015 
0.854 ±0.033 
1.048 ±0.026 
0.742 ± 0.046 
0.0610 ±0.0019 
0.03046 ± 0.00051 



9.9 



+3.8 



151.0 
0.778 
0.798 
4.525 

0.829 
1.033 
0.753 
0.0619 
0.03002 
15. 



±1.9 
±0.029 
±0.014 
±0.015 
±0.026 
± 0.023 
± 0.046 
±0.0018 
± 0.00037 

1+2.6 
^-2.3 



151.2±3. 
0.782 ±0. 
0.799 ±0. 
4.526 ±0. 
0.832 ±0, 
1.034 ±0. 
0.752 ±0. 
0.0618 ±0. 

0.03006 ± 0. 
14.812 





047 
019 
017 
038 
029 
048 
0021 

00059 

6 

4 



152.8 ± 
0.807 ± 
0.808 ± 
4.530 ± 
0.849 ± 
1.045 ± 
0.744 ± 
0.0612 ± 
0.03038 ± 
9.3 



2.2 
0.035 
0.015 
0.015 
0.030 
0.025 
0.046 
0.0019 
0.00044 

2.6 
2.7 



data would further improve measurements of the physical proper- 
ties of this planetary system. 
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Table 10. Final physical properties of the WASP-2 system plus a comparison to literature measurements. 



This work (final) 



Charbonneau Torres et al. 
et al. (2007) (2008) 



Southworth (2009) 



Daemgen et al. Triaud et al. 
(2009) (2010) 



MaCMq) 
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0.813 ±0.032 


Ra{Rq) 


0.807 ± 0.019 ±0.011 


0.81 ±0.04 


log 9a (cgs) 


H^.DDyj It U.Ui / It U.UUO 




PA (P©) 


1.531 ±0.067 




Mb (Mjup) 


0.847 ± 0.038 ±0.024 




(Rjup) 


1.044 ± 0.029 ±0.015 


1.038 ±0.050 


3b (ms-1) 


19.30 ±0.82 




Pb ( Pjup) 


0.745 ± 0.048 ±0.010 




TJ^ (K) 


1280 ±21 







0.0612 ± 0.0021 ±0.0009 




a(AU) 
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Age (Gyr) 


11 S+8-1 +3-3 
lJ^-» _4.4 -2.5 





0.89t«;- 
0.840in^^ 

7+0.035 



4.537 
1.45 
0.915 
1.071 



0.046 
+0.19 
0.11 
+0.090 
0.093 
+0.080 
0.083 



19.4t^^ 

74+0-22 
^•'^-0.16 

13041^4 

0592+° °°"'® 
+0.00130 
0.00154 
+8.4 
5.6 



0.03138 
5.6 
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0.03120 ±0.00130 ±0.00030 
unconstrained 



0.89 ±0.12 



1.49 ±0.15 



Q 04+0.11 
^•"^-0.12 



0.843 ±0.063 0.825+q;o4o 



1.491 



+0.088 



0.914 ±0.092 0.866lg:o8| 



1.117 ±0.082 1.077; 
19.0± 1.6 
0.70 ±0.19 
1300 ±54 
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-0.0015 
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